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Abstract 

We investigate the structure of equations of motion and lagrangian constraints in 
a general theory of massive spin 2 field interacting with external gravity. We demon- 
strate how consistency with the flat limit can be achieved in a number of specific 
spacetimes. One such example is an arbitrary static spacetime though equations 
of motion in this case may lack causal properties. Another example is provided 
by external gravity fulfilling vacuum Einstein equations with arbitrary cosmological 
constant. In the latter case there exists one-parameter family of theories describing 
causal propagation of the correct number of degrees of freedom for the massive spin 
2 field in arbitrary dimension. For a specific value of the parameter a gauge invari- 
ance with a vector parameter appears, this value is interpreted as massless limit 
of the theory. Another specific value of the parameter produces gauge invariance 
with a scalar parameter and this cannot be interpreted as a consistent massive or 
massless theory. 
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Problems of consistent equations of motion for interacting higher spin fields deserve 
studying due to many reasons. First of all, string theory includes an infinite tower of 
massive excitations with all possible spins and thus should allow some consistent effective 
description of arbitrary spin fields interaction. Second, composite resonance particles with 
higher spins do exist and one should be able to describe their interaction (for example, 
with external electromagnetic and gravitational fields) in terms of some effective local 
field theory. At last, investigation of higher spin fields is interesting on its own from the 
general point of view. It would be surprising if nature admits description of free fields 
with arbitrary spins but stops, say, at spin 1 in case of interacting massive fields. Even if 
it is really the case, one should try to understand why this is the way the nature works. 

This note is devoted to investigation of the massive spin 2 field interacting with exter- 
nal gravity which represents one of the simplest higher spin models. It has been studied in 
numerous papers but careful and general analysis of the consistency and causality 

of the theory in arbitrary curved spacetime was still absent. 

There are at least two ways the interaction may spoil the consistency of a higher spin 
theory. Firstly, interaction may change the number of dynamical degrees of freedom. For 
example, a massive field with spin s in D = 4 Minkowski spacetime is described by a rank 
s symmetric traceless transverse tensor satisfying the mass shell condition: 

(d 2 - m 2 )0 M1 ... Ms = 0, dy^...^ = 0, 0V,..^_ 2 = 0. (1) 

To reproduce all these equations from a single lagrangian one needs to introduce auxiliary 
fields Xm...n a -ai Xm—nssi ■ ■ ■ > X HH 0- These symmetric traceless fields vanish on shell 
but their presence in the theory provides lagrangian description of the conditions (|T|). In 
higher dimensional spacetimes there appear fields of more complex tensor structure but 
general situation remains the same, i.e. lagrangian description always requires presence 
of unphysical auxiliary degrees of freedom. 

Namely these auxiliary fields create problems when one tries to turn on interaction in 
the theory. Arbitrary interaction makes the auxiliary fields dynamical thus increasing the 
number of degrees of freedom. Usually these extra degrees of freedom are ghostlike and 
should be considered as pathological. Requirement of absence of these extra dynamical 
degrees of freedom imposes severe restrictions on the possible interaction. 

To construct a consistent massive field theory one often starts with a corresponding 
consistent massless theory which should be invariant with respect to gauge symmetry 
and then breaks this symmetry by introducing mass terms into lagrangian. In this case 
invariance of the kinetic part of lagrangian guarantees the correct number of degrees of 
freedom in both massless and massive theories. For example, spin 2 field possesses a gauge 
invariant massless lagrangian only if external gravity satisfies vacuum Einstein equations 
and so it is usually believed that the massive spin 2 field can consistently propagate also 
only in Einstein spacetimes. 

In this paper which is a sequel to [[H]] we show that this belief is not true and that 
in massive case there exists a number of possibilities of providing the correct number 
of degrees of freedom. Of course, invariance of kinetic term does always provide the 
correct number of degrees of freedom but this is not the only possibility. In particular, 
we describe below an example when external spacetime does not fulfill Einstein equations 
but the theory is consistent with the flat spacetime limit. 

2 For analysis of the corresponding massless models see e.g. |j| ||, |ij 
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The other problem that may arise in higher spin fields theories is connected with 
possible violation of causal properties. This problem was first noted in the theory of spin 



3/2 field in external fields [|T3| (see also the review \TA\ and a recent discussion in ||15|| ) 

In general, when one has a system of differential equations for a set of fields <fi B (to be 
specific, let us say about second order equations) 

MAB^d^A 3 + • • • = 0, fi, v = 0, . . . , D - 1 (2) 

the following definitions are used. A characteristic matrix is the matrix function of D 
arguments n p built out of the coefficients at the second derivatives in the equations: 
MAB{ n ) — MAB^ Un [i n v A characteristic equation is det M^^n) = 0. A characteristic 
surface is the surface S(x) = const where d^S(x) = n M . 

If for any nj (i = 1, . . . , D—l) all solutions of the characteristic equation no(rii) are real 
then the system of differential equations is called hyperbolic and describes propagation of 
some wave processes. The hyperbolic system is called causal if there is no timelike vectors 
among solutions n M of the characteristic equations. Such a system describes propagation 
with a velocity not exceeding the speed of light. If there exist timelike solutions for 
then the corresponding characteristic surfaces are spacelike and violate causality. 

Turning on interaction in theories of higher spin fields in general changes the character- 
istic matrix and there appears possibility of superluminal propagation. Such a situation 
also should be considered as pathological. 

Both these problems arise in the theory of massive spin 2 field coupled to external 
gravitational field. To provide consistency of the interaction we should conserve the same 
number of physical degrees of freedom and constraints that the theory possesses in flat 
spacetime. To find the complete set of constraints we will use the general lagrangian 



scheme |L6[ which is equivalent to the Dirac-Bergmann procedure in hamiltonian formal- 
ism but for our purposes is simpler. In the case of second class constraints (which is 
relevant for massive higher spin fields) it consists in the following steps. If in a theory of 
some set of fields (p A {x), A = 1, . . . , N the original lagrangian equations of motion define 
only r < N of the second time derivatives ( "accelerations" ) 4> A then one can build N — r 
primary constraints, i.e. linear combinations of the equations of motion that does not con- 
tain accelerations. Requirement of conservation in time of the primary constraints either 
define some of the missing accelerations or lead to new (secondary) constraints. Then one 
demands conservation of the secondary constraints and so on, until all the accelerations 
are defined and the procedure closes up. 

Before considering the theory in external gravitational field we analyze the structure 
of equations of motion in Minkowski spacetime. The purpose of this analysis is twofold. 
First, we illustrate the general scheme of calculating the constraints within covariant 
lagrangian framework. In addition, building a consistent and causal theory in curved 
spacetime we use these flat constraints as a reference point. 

Free spin 2 field is known to be described by the Fierz-Pauli action [|ITJ (we consider 
arbitrary spacetime dimension): 

$ = I d D x{- A d^Hd»H - l -d,H vp d»H»P - l -d»H, v d»H + ^H up d?H^ 

i i 
m m 9 

T z 



-—H flu H^ + —H^. (3) 
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Here the role of auxiliary field is played by the trace H = rf^H^. The equations of 
motion 

= d 2 H^- V ^d 2 H + d^H + V ^d a dPH aP -d a d^H\-d a d u H% 

- m 2 H^ + m 2 Hrj^ = (4) 

contain D primary constraints (expressions without second time derivatives H^ v ): 

E 00 = AH«-d i d j H ij -m 2 H ii = <p$ ) KiO (5) 
E 0i = AH 0i + diH kk - d k H ki - did k H ok - m 2 H 0i = ipf 1 « 0. (6) 

The remaining equations of motion Eij = allow to define the accelerations i/jj in terms 
of H^y and H^ u . The accelerations H 00 , H 0i cannot be expressed from the equations 
directly. 

Conditions of conservation of the primary constraints in time E 0fl « lead to D 
secondary constraints. On-shell they are equivalent to 

ipW = d^E^ = m 2 d u H - m^H^ w (7) 

Conservation of (p\ 2 ^ defines D — 1 accelerations H^i and conservation of (p^ gives 
another one constraint. It is convenient to choose it in the covariant form by adding 
suitable terms proportional to the equations of motion: 

(3) = Qnff> E + J^^e = Hm ±2^1 ~ o (8) 
y i» D _ 2 / fw D -2 v ' 

Conservation of gives one more constraint on initial values 

y? (4) = -J?oo + H kk = H « (9) 

and from the conservation of this last constraint the acceleration ifoo is defined. 

Altogether there are 2D + 2 constraints on the initial values of H^ v and H^ v . The 
lagrangian theory is equivalent to the system of the equations 

(a 2 -m 2 )#^ = 0, 3^ = 0, H% = 0. (10) 

and describes traceless and transverse symmetric tensor field of the second rank. 

Obviously, the equations of motion ( |T0| ) are causal because the characteristic equation 

detM(n) = (n 2 ) D ^/ 2 (11) 

has 2 multiply degenerate roots 

-rag + n? = 0, n = ±Vn?. (12) 

which correspond to real null solutions for n^. Note that analysis of causality is possible 
only after calculation of all the constraints. Original lagrangian equations of motion @ 
have degenerate characteristic matrix det M(n) = and do not allow to define propagation 
cones of the field H^ u . 
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In the massless limit m 2 = the structure of the theory (|3]) changes. Instead of 
the secondary constraints (0) conservation of the primary constraints lead to identities 
d^E^ = which mean that the theory becomes gauge invariant with respect to the 
local transformations 5H^ U = d^ v + d u ^^. Such a theory represents the quadratic part 
of the Einstein-Hilbert action for gravitational field and the gauge invariance is a linear 
counterpart of the general coordinate invariance. 

Now if we want to construct a theory of massive spin 2 field on a curved manifold 
first of all we should provide the same number of propagating degrees of freedom as in 
the flat case. It means that new equations of motion E^ u should lead to exactly 2D + 2 
constraints and in the flat spacetime limit these constraints should reduce to their flat 
counterparts. 

Generalizing (||) to curved spacetime we should substitute all derivatives by the co- 
variant ones and also we can add non-minimal terms containing curvature tensor with 
some dimensionless coefficients in front of them. As a result, the most general action for 
massive spin 2 field in curved spacetime quadratic in derivatives and consistent with the 
flat limit should have the form ||T|: 

8 = J d D xV^G~\^V fl HV»H - jV IM H up V»H 1/ o - ^WH^VH + ^V^H^H^ 

, fll p tt zraf3 I ° 2 D fJ-2 | °3 -Qiiavfi tt tt , a 4 pa/3 jj tt a , a 5 tt tt 

+ —tiii a pti + + -tlfiu-tlap + —tt n aa np + —K U a/3 H 

-^H^ + ^H 2 ] (13) 
where ax, . . . as are so far arbitrary dimensionless coefficients, R^uXk = d\T£ K — . . ., R^ u = 
Equations of motion 

+ 2a\RH iiV + 2a2G llv RH + 2aj,R^ a v H a p + a^R^H^ + a^R^H^ 

+ asR^H + a b G^R aP H aj3 - m 2 H IIU + m 2 HG llu « (14) 

contain second time derivatives of in the following way: 

-£%) — [vr — LrQoLr Gr + Lr 00 Lr (jr ) V Q V O-H-mn + U [ V o) , 

E 0l = (-G 0i G 00 G m " + G' 0l G' 0m G "-G' 0m ( 5r)VoVo// m „ + 0(Vo), 

Eij = (G 00 5™5j - GijG 00 G mn + GijG 0m G 0n )VoV H mn + O(V ). (15) 

So we see that accelerations H 00 and H 0i again (as in the flat case) do not enter the 
equations of motion while accelerations Hij can be expressed through H^, and their 
spatial derivatives. 

There are D linear combinations of the equations of motion which do not contain 
second time derivatives and so represent primary constraints of the theory: 

<p« = E\ = G 00 E , + G 0j E jfM (16) 

Now one should calculate time derivatives of these constraints and define secondary ones. 
In order to do this in a covariant form we can add to the time derivative of (pft any linear 
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combination of equations of motion and primary constraints. So we choose the secondary 
constraints in the following way: 



(2aii2 - m 2 )V»H^ + (2a 2 i? + m 2 ) V V H + 2a 3 ^ a /V M ^ + a 4 R^ a V^H a 
+ (04 - 2)R a u V»H afi + a 5 R a ^V v H afl + (a 5 + l)R a v V a H 

+ (2a, + ^)H au V a R + (2a 2 + °^)HV V R 



a/3 



(2a 3 + a 5 + l)V,E a ^ + (a 4 - 2a 3 - 2)V a i? 



a p/3 



(17) 



At the next step conservation of these D secondary constraints should lead to one new 
constraint and to expressions for D—l accelerations H Qi . This means that the constraints 
( PD should contain the first time derivatives H 0fl through the matrix with the rank D—l: 



(2) 

</>o 
,( 2 ) 



rank = rank 
In the flat spacetime we had the matrix 



$ v 



a^oo + DjHoj + ... 



A 




a 


Dj 



D 



(18) 
(19) 





m 2 5i 



(20) 



In the curved case the explicit form of this matrix elements in the constraints ( |T7| ) is: 



.4 
B j 

Ci 



RG m (2 ai + 2a 2 ) + R m (a 4 + a 5 ) + iT G uu (a 4 + a 5 - 1) 
m 2 G 0j + RG 0j {2 ai + 4a 2 ) + 2a 3 R 0j ° + R j G 00 {a 4 - 2) 
+ R 0j (a 4 + 2a 5 ) + R° G 0j (a 4 + 2a 5 ) 
R°iG m (a A + a^-l) 



m 



+ «5 

'G m 5{ + 2a 1 RG 00 5i + 2a 3 iW + a 4J R 00 ^' 



+ (a 4 - 2)R j iG 00 + (a 4 + 2a 5 )R°G 0j 



(21) 



At this stage the restrictions that consistency imposes on the type of interaction reduce 
to the requirements that the above matrix elements give 



det $ = 0, det ^ 



(221 



When the gravitational background is arbitrary it is not clear how to fulfill this condition 
by choosing some specific values of non- minimal couplings a±, . . . a 5 . For example, require- 
ment of vanishing of the elements A and Q ([H]) would lead to contradictory equations 
a 4 + a 5 = 0, a 4 + — 1 = 0. 

But the consistency conditions ( $2\) can be fulfilled in a number of specific gravitational 
background. Namely, any spacetime which in some coordinates has 



R°i = 



(23) 



5 



provides such an example. In such a spacetime R 00 = R°oG 00 and choosing coefficients 
&i + 0-2 — 0, 2a 4 + 2a5 = 1 we have the first column of the matrix $ vanishing and so the 
conditions (^2|) fulfilled. 

As a first example where ( p3|) holds let us consider an arbitrary static spacetime, i.e. a 
spacetime having a timelike Killing vector and invariant with respect to the time reversal 



x 



-x°. In such a spacetime one can always find coordinates where 



0, G 0i = 0. 
The matrix elements (BID in this case become 



(24) 



A = J RG 00 (2ai + 2a 2 ) + ^ 00 (2a4 + 2a 5 -l), B j = 0, Q = 0, 
Dj = ( - m 2 G 00 + 2 ai RG 00 + a A R m )5l + (a 4 - 2)i^G 00 + 2a 3 R 0j i ° 



and (|22 ) lead to the following conditions: 

2ai + 2a 2 = 0, 2a 4 + 2a 5 - 1 = 0, det iy ^ 



(25) 



(26) 



The last inequality may be violated in strong gravitational field and as we comment below 
this fact may lead to causal problems. 

Suppose that all the conditions fl26|) are fulfilled. For simplicity we also choose a 3 = 0. 
Then we have the equations of the form fllil) with the coefficients 
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~, «3 = 0, a 4 = i - £ 2 , a 5 = 6 



(27) 



where £i, 6 are two arbitrary coupling parameters. 
One of the secondary constraints 



(2) 



V a E, 



aO 



G i3 [m 2 - i x R) + (1 + &)G ij R% + &R ij 



VoHij 



GV(Z l R-m')-[- + &)GVR u 



6if 



6 )-^OiV i i?°o 



(28) 



does not contain velocities H 00 , H 0i and so its conservation leads to a new constraint 
y?( 3 ) »i VoV a -E' a o- After exclusion from this expression the accelerations f/y we get this 
constraint as the following combination of the equations of motion: 



(3) 



D-2 



m 2 G 00 + (6 - fi)i2G o + R 



00 



D — 1 
D-2 



m + 



26-26(^-1) 



D-2 



D — 1 
D-2 



m 2 R°o - QRijR 1 - 



+ gRR + 



D-2 



2^6 )RR°o + 6(6 + 1)^00^°° #00 + . . .(29) 



We did not write down the explicit form of this constraint because everything we 
should know about it is the way it contains the component Hqq. Namely, contains 
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neither the acceleration Hqq nor the velocity ifoo- It means that its conservation in time 
leads to another new constraints 

V? {4) « V ^ (3) (30) 

and hence the total number of constraints is the same as in the flat spacetime provided 
that the expression in the braces in front of Hqq in tp^ does not vanish. 

Vanishing of this expression in braces as well as violation of the inequality (|26| ) leads to 
local changing of the number of degrees of freedom and this fact is known to be related with 
acausal behavior in higher spin theories in external fields |TR 117). In general, causality 



breaks in those cases when there are points in spacetime in which it is impossible to define 
all the accelerations from the conservation of constraints. 

In our case it means that causality will hold everywhere only if det D l j ^ and the 
expression in braces in tp^ also does not vanish. Obviously, in general case there are 
values of R^ that violate these requirements. 

It is instructive to consider in more detail the Reissner-Nordstrom solution in D = 4 
simple example of non-trivial static spacetime 



tf = -(l-^ + g)tf + - g-y + W (31) 



In this case causality problems are absent when the expressions 



(3) _ {£_l m 4 +m2 ( 26 + £_i)5! +6(1 _ 26 )^} ffoo + ... ( 32) 



do not vanish. Far enough from the horizon where all the terms containing r in ( |32D are 
too small and so propagation is causal in this region. Causal problems may develop only 
for small values of r. This might be excluded if all terms in fl32|) were positive, that is if 



l+2£ 2 <0, 2 £ 2 + ^-|>0, k(l-2&)>0 (33) 

but these three conditions are contradictory. 

It means that for any value of the coupling parameter expressions in fl32|) vanish for 
some values of the coordinate r and the massive spin 2 field propagate causally only in 
the regions near infinity but close to the horizon causality is lost. Of course, this example 
does not mean that there cannot exist other spacetimes where causality might be achieved 
everywhere for some special values of coupling parameters. 

Another possible way to fulfill the consistency requirements ( p2|) is to consider space- 
times representing solutions of vacuum Einstein equations with arbitrary cosmological 
constant: 

R^u — j^G^R . (34) 
In this case the coefficients a^, 05 in the lagrangian (|13"D are absent and the matrix $ 
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takes the form: 



$ v = 



i?G 00 (2a 1 + 2a 2 - jj) 


i2G° J '(2oi + 4a 2 ) + 2a 3 R 0j ° + m 2 G 0j 





2a 3 R°h° + RG 00 Sl(2a 1 - %) ~ m 2 G 00 5j 



(35) 



The simplest way to make the rank of this matrix to be equal to D — 1 is provided by the 
following choice of the coefficients: 



2ai + 2a 2 - = 0, a 3 = 0, 



0. 



(36) 



As a result, we have one-parameter family of theories 



ai= D' ° 2= 2D ' 



03 = 0, a 4 = 0, as = 



i2 



'/XI/ 



1 



(37) 



with £ an arbitrary real number. 

The action in this case takes the form 



(38) 



-H' 



and the corresponding equations of motion are 

(39) 



2£ 1 — 2^ 
+ j^RH flu H — — RHG^ 



m 2 H^ + m 2 HG fMU = 



The secondary constraints built out of them are 

= W a E aix = (V,H - V a H, a ) (m 2 + 
and the matrix $ looks like 



2(1-0 



D 



R 



2(1-0 














(40) 



(41) 



(2) 

Just like in the flat case, in this theory the conditions (p\ ~ define the accelerations 

(2) 

and the condition ^' «0 after excluding if i gives a new constraint, i.e. the acceleration 
.ffoo is n °t defined at this stage. 

To define the new constraint in a covariant form we use the following linear combination 
of (fiff, equations of motion, primary and secondary constraints: 



<P 



(3) 



H 



D-2 



G^E ul/ + VTB F 4 
R + m 



2(1-0 



D-2~ ^ 

1 /2(l-0 



D 



"" I HI) 2) "" 
2 \(D + 2(j{l-D) r> , 
D 



RG^E^ = 

R + m 2 (D-l)\ w0. (42) 
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This gives tracelessness condition for the field provided that parameters of the theory 
fulfill the conditions: 

?S^« fl + m V0, D + 2 f- D) R + m '(D-l)^0 (43) 

Requirement of conservation of (p^ leads to one more constraint 

£( 3 ) ~ H <p {A) = H « 0. (44) 

The last acceleration ifoo is expressed from the condition <^( 4 ) « 0. 

Using the constraints for simplifying the equations of motion we see that the original 
equations are equivalent to the following system: 

V 2 H^ V + 2R a / v H aP + 2(C ~ l) RH^ v - m 2 H flu = 0, 

H% = 0, H^ = 0, V^ = 0, (45) 

G 00 V V l H\ - G 0i V V^°, - G^V.Votf , - G ij ViVjH Q v - 2R a0 ^H af3 

- 2 -^lRH\ + m 2 H\ = 0. 
The last expression represents D primary constraints. 



For any values of £ (except two degenerate values excluded by ((43])) the theory de- 
scribes the same number of degrees of freedom as in the flat case - the symmetric, covari- 
antly transverse and traceless tensor. D primary constraints guarantees conservation of 
the transversality conditions in time. 

Let us now consider the causal properties of the theory. Again, if we tried to use the 
equations of motion in the original lagrangian form (|39"D then the characteristic matrix 



M^ A » = <V) {AK V - G, u G XK n 2 + G XK n,n v + G^nV - 5^n^ - 5^n v (46) 

would be degenerate. This fact can be seen from the relation 

n»M^ XK (n) = (47) 

which means that any symmetric tensor of the form n^t u ^ (with t v an arbitrary vector) 
represents a "null vector" for the matrix M(n) and therefore det M — 0. 

After having used the constraints we obtain the equations of motion written in the 
form fl4"5| ) and the characteristic matrix becomes non-degenerate: 

M^ XK (n) = <V A V, n 2 = G a ^n a n p . (48) 

The characteristic cones remains the same as in the flat case. At any point x Q we can 
choose locally G Q/3 (xo) = r] af3 and then 



n 2 



J-0 



-n 2 + nf (49) 



Just like in the flat case the equations are hyperbolic and causal. 

Now let us discuss the massless limit of the theory under consideration. There are sev- 
eral points of view on the definition of masslessness in a curved spacetime of an arbitrary 
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dimension. We guess that the most physically accepted definition is the one referring to 
appearance of a gauge invariance for some specific values of the theory parameters (see 



e.g. |fT9], [20] for a recent discussion). 

In our case it means that the real mass parameter M for the field H^ v in an Einstein 
spacetime is defined as 

M 2 = m 2 + 2{ - l D (50) 

When M 2 = instead of D secondary constraints ip^ we have D identities for the 
equations of motion V^E^ = and the theory acquires gauge invariance 5H^ U = + 
V,/^. This explains the meaning of the first condition in (|43|), it just tells us that the 
theory is massive. 

In fact, two parameters m 2 and £ enter the action (^) in a single combination M 2 
([50]). Since scalar curvature is constant in Einstein spacetime there is no way to distinguish 
between the corresponding terms ~ £RHH, ~ m 2 HH (with arbitrary £, m) in the action. 
The difference between the two will appear only if we consider Weyl rescaling of the metric. 
Note that the "massless" theory with M 2 = is not Weyl invariant. In the case of dS/AdS 
spacetimes the difference between masslessness, conformal and gauge invariance and null 
cone propagation was discussed in detail in ET|. In our case the theory obviously cannot 



possess Weyl invariance. 

The second inequality Q35D is more mysterious. If it fails to hold, i.e. if M 2 = M 2 
D ~%y\R then instead of the constraint p^ the scalar identity 



D(D-l)- 



R_ 
D(D~~ 1) 



VT£ F + ———G^E^ = (51) 



with the corresponding gauge invariance 

R 

5H^ = V M V,e + D r D _ x) G^e (52) 

arise. 

Appearance of this gauge invariance with a scalar parameter was first found for the 



massive spin 2 in spacetime of constant curvature in |21] and was further investigated 0, |J] 
in spacetimes with positive cosmo logical constant. Our analysis shows that this gauge 
invariance is a feature of more general spin 2 theories in arbitrary Einstein spacetimes. In 
this case we can simplify the equations of motion using the secondary constraints (pE0|): 

V 2 H^ - V^J U H + 2R^/H a/3 + ^ RH^ ~ ^ RG^H = 0. (53) 

After imposing the gauge condition^ H = one can see that these equations describe 
causal propagation of the field but the number of propagating degrees of freedom 
corresponds to neither massive nor massless spin 2 free field. It was argued in 0, || that 
appearance of the gauge invariance ( |5*2]) leads to such pathological properties as violation 
of the classical Hamiltonian positiveness and negative norm states in the quantum version 



3 It does not fix (|5^) completely and the residual symmetry with the prameter obeying (V 2 + -p^y)e = 
remains. 
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of the theory. One should expect similar problems in the general spin 2 theory in arbitrary 
Einstein spacetime described in this paper. 

We demonstrated that correct number of degrees of freedom in the massive spin 2 
theory (algebraic consistency) can be achieved in a large class of curved spacetimes which 
include as particular cases arbitrary static spacetimes and vacuum Einstein spacetimes. 
An analysis of the constraints structure shows that in case of a static spacetime there 
esxists a potential source of acausal behavior. However, as we see on the example of 
Reissner-Norsdtrom spacetime causal propagation is possible in the regions where gravita- 
tional field is weak enough. In Einstein spacetimes spin 2 massive field can be consistently 
described by a one-parameter family of theories ( |38| ) . For any value of the parameter sat- 
isfying (P^|) the corresponding equations describe the correct number of degrees of freedom 
which propagate causally. 



It is interesting to compare our approach with the paper [22[] which is devoted to 



investigation of consistency of higher rank spin-tensor fields in curved spacetime from a 
different point of view. The authors of [22f considered equations for the fields carrying 
irreducible representations of Euclidian version of the four dimensional Lorentz group 
SO (4) and analyzed when irreducibility of these representations is preserved in curved 
space. In particular, they showed that symmetric second rank tensor equations are con- 
sistent in this sense in Einstein spaces. However, such an analysis (sufficient for the proof 
of index theorems in P%| ) is not enough when one tries to build a consistent theory for 
a physical field on a curved manifold starting from an irreducible representation of the 
Poincare group in flat spacetime. Preservation of the correct number of degrees of freedom 
in such a theory is a requirement independent from the algebraic consistency considered 
in \22\. Einstein spacetimes provide an example when both these conditions are fulfilled 
but as we saw in our analysis correct number of degrees of freedom can be preserved in 
much wider class of spacetimes. Besides, in physical theories for interacting higher spin 
fields we face a new problem of causality which should also be studied independently. 
In general there can be theories with correct number of degrees of freedom but acausal, 
and we really see such examples in case of spin 2 field in external gravity. It is worth to 
note that the authors of emphasized that their analysis has no direct relation to the 
problem of consistent propagation of higher spin physical fields and that they did not set 
this problem at all. 

In case of Einstein spacetimes our lagrangian for the spin 2 field in curved spacetime 
is the most general known so far, in all previous works only the theories with specific 
values of the parameter £ were considered [f|, H . Two degenerate values of the parameter 
£ describe the theories with different degrees of freedom. One of this degenerate values 
corresponds to massless spin 2 field in an Einstein spacetime, another one describes neither 
massive nor massless spin 2 field. 

The next natural step would consist in building a theory describing dynamics of both 
gravity and massive spin 2 field. In such a theory in addition to dynamical equations for 
the massive spin 2 field one would have dynamical equations for gravity with the energy- 
momentum tensor constructed out of spin 2 field components. The analysis of consistency 
then changes and one needs to have correct number of constraints and causality for both 
fields interacting with each other |]]. 

The only known consistent system of a higher spin field interacting with dynamical 
gravity is the theory of massless helicity 3/2 field, i.e. supergravity (see also the book 
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||24|| ). In that case consistency with dynamical gravity requires four-fermion interaction. 
If a consistent description of spin 2 field interacting with dynamical gravity exists it may 
also require some non-trivial modification of the lagrangian. At least, it is known that 
lagrangians quadratic in spin 2 field do not provide such a consistency |J. A possible 
way of consistent description of the spin 2 field on arbitrary gravitational background 
was recently proposed in |K|. This was achieved by means of representation of the la- 



grangian in the form of infinite series in curvature and imposing the consistency condition 
perturbatively in each order (earlier similar construction was investigated for symmetric 
Einstein spacetime in [El). 

Further generalizations of our analysis may include theories of massive spins s > 3 
fields (which would require more complex structure of auxiliary fields) and interaction 
with other background fields, e.g. with scalar dilaton and antisymmetric tensor that are 
relevant in string theory. 
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